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O ; ABSTRACT 

We present simultaneous optical and X-ray data from Swift for a sample of radio-loud flat 
spectrum quasars selected from the Einstein Medium Sensitivity Survey (EMSS). We present 
also a complete analysis of Swift and INTEGRAL data on 4 blazars recently discussed as 
. possibly challenging the trends of the hypothesised "blazar spectral sequence". The SEDs 

Qh' of all these objects are modelled in terms of a general theoretical scheme, applicable to all 

blazars, leading to an estimate of the jets' physical parameters. Our results show that, in 
j!^ i the case of the EMSS broad line blazars. X-ray selection does not lead to find sources with 

^ ' synchrotron peaks in the UV/X-ray range, as was the case for X-ray selected BL Lacs. Instead, 

Ci ', for a wide range of radio powers all the sources with broad emission lines show similar SEDs, 

with synchrotron components peaking below the optical/UV range. The SED models suggest 
that the associated IC emission should peak below the GeV range, but could be detectable 
in some cases by the Fermi Gamma-Ray Space Telescope. Of the remaining 4 "anomalous" 
1^ ■ blazars, two highly luminous sources with broad lines, claimed to possibly emit synchrotron 

■r^-j- I X-rays, are shown to be better described with IC models for their X-ray emission. For one 

source with weak emission lines (a BL Lac object) a synchrotron peak in the soft X-ray range 
' is confirmed, while for the fourth source, exhibiting narrow emission lines typical of NLSy Is, 

no evidence of X-ray emission from a relativistic jet is found. We reexamine the standing and 
interpretation of the original "blazar spectral sequence" and suggest that the photon ambient, 
' in which the particle acceleration and emission occur, is likely the main factor determining 

, the shape of the blazar SED. A connection between SED shape and jet power/luminosity 

can however result through the link between the mass and accretion rate of the central black 
hole and the radiative efficiency of the resulting accretion flow, thus involving at least two 
, parameters. 
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1 INTRODUCTION 

The Spectral Energy Distributions (SED) of blazars from radio to 
gamma-rays exhibit remarkable properties, in particular a "univer- 
sal" structure consisting of two broad humps. A systematic com- 
pilation of multiwavelength data for three complete samples, two 
selected in the radio band (2Jy FSRQ and 1 Jy BL Lacs) and one 
in X-rays (Einstein Slew Survey BL Lacs) was performed by Fos- 
sati et al. (1998). Note that, at the time, the results from the Comp- 
ton Gamma-Ray Observatory (CGRO) were still under study, while 
very few objects had been detected in the TeV band, thus the data 
at high energies were rather incomplete. 

The derived SEDs were averaged in fixed radio luminosity 
bins, irrespective of sample membership or classification (BL Lac 
vs. FSRQ). The result was the well known "Spectral Sequence" (a 
series of 5 average SEDs, "sequence" in the following) showing 
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the peaks of the two spectral components shifting systematically to 
higher frequencies with decreasing luminosity. For brevity we will 
refer to SEDs with synchrotron peak below or above the optical 
range as "red" or "blue" SEDs respectively. 

It has to be stressed that the derived average SEDs, by them- 
selves, only represent the spectral properties of the considered sam- 
ples and do not contain further implications; however the fact that, 
despite the mixing of different samples and classifications, the av- 
erage SEDs show systematic trends, possibly driven by the radio 
power suggests an underlying regularity of behaviour of jets in 
very different objects. Understanding this behavior could lead to 
a more profound comprehension of the jets' radiation mechanisms 
and physics. 

Ghisellini et al. (1998) modelled a large number of individual 
sources with sufficient data at high energies with a simple general 
scheme: a single emission region, moving with bulk Lorentz fac- 
tor r ~ 10, contains relativistic electrons emitting via synchrotron 
radiation the lower energy spectral component; the same electrons 
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up scatter synchrotron photons as well as photons from an exter- 
nal radiation field producing the high energy spectral component. 
The electron spectrum is computed assimring a broken power law 
injection spectrum and taking into account radiative energy losses. 

With these minimal assumptions the "sequence" of spectral 
properties translates into a parameter sequence, in which the en- 
ergy of the electrons radiating at the SED peaks is higher for lower 
values of the (comoving) total magnetic field plus photon energy 
densities, due to the reduced energy losses. This interpretation of 
the "sequence" concept and following developments (Ghisellini, 
Celotti & Costamante 2002; Costamante & Ghiselhni 2002) led to 
successful predictions of the best candidates for TeV detection (see 
Wagner 2008 for a recent review). 

On the other hand, the phenomenological "sequence" is likely 
affected by observational biases (Maraschi and Tavecchio 2001 
(ASP CS Vol 227)). Significant efforts were made to construct new 
samples of blazars with different thresholds and selection criteria 
{e.g., Laurent-Muehleisen et al. 1999, Caccianiga & Marcha 2004, 
Massaro et al. 2005, Landt et al. 2006, Giommi et al. 2007a, Tur- 
riziani et al. 2007, Healey et al. 2008, Massaro et al. 2008) and 
some discrepant objects were found (see Padovani 2007 for a re- 
view on the challenges to the blazar sequence). In particular Cac- 
cianiga & Marcha (2004) evidenced the existence of sources with 
relatively low radio power, but X-ray to radio ratios similar to those 
ofFSRQs. 

In a somewhat different line of approach the group lead by 
Valtaoja developped a technique to estimate the beaming parame- 
ters (Doppler factors Lorentz factors and viewing angles) of blazar 
jets from the radio properties alone (Valtaoja et al. 1999, Lahteen- 
maki, Valtaoja & Wiik 1999, Lahteenmaki & Valtaoja 1999). Bas- 
ing on extensive monitoring data at high radio frequencies (22 and 
37 GHz) these authors compute "apparent" brightness temperatures 
using light travel times inferred from variability as sizes for the 
emission region. Comparing the results with the theoretical maxi- 
mvun brightness temperature allowed by the SSC process they de- 
rive the Doppler factor of the emission region and further use the 
apparent superluminal speeds to determine the bulk Lorentz factor 
and the viewing angle. This procedure is correct, however one has 
to keep in mind that the results obtained refer to the radio emitting 
regions of the jet, which do not coincide with the regions emit- 
ting the high energy radiation (e.g. Ghisellini et al. 1998; Sikora 
2001). Moreover in the analysis of the SEDs of a large number of 
BL Lac objects Nieppola et al. (2006) fit the data in the full (radio 
to gamma-ray) frequency range available with a single paraboUc 
curve, thus their determination of the "peak frequency" of the SED 
cannot be compared with our approach. In Nieppola et al. (2008) 
however, the parabolic fits are not forced to include the X-ray and 
higher energy data, yielding peak frequencies that can be compared 
with ours. We will discuss their results in the Discussion Section. 

One of the problems in the "sequence" assembly was the lack 
of an X-ray selected sample of flat-spectrum radio quasars. Wolter 
& Celotti (2001) selected and discussed a sample of FSRQ within 
the X-ray selected sample of radio Loud AGN from the Einstein 
Medium Sensitivity Survey (EMSS). They found that the average 
broad-band spectral indices of these X-ray selected FSRQs were 
consistent with those of the radio selected FSRQs within the "se- 
quence". However, a study of their SEDs was not possible due to 
insufficient spectral data. 

Here, we present Swift observations of the 10 X-ray bright- 
est objects in the Wolter & Celotti sample as well as a full analy- 
sis of Swift data for 4 blazars recently claimed to be at odds with 
the sequence scheme. We model the SEDs of the 14 objects fol- 



lowing Ghisellini, Celotti & Costamante (2002), deriving estimates 
of important physical quantities for their jets. Finally, in the light 
of these new results, we discuss the present standing and interpre- 
tation of the blazar "spectral sequence". Preliminary results were 
given in Maraschi et al. (2008). The scheme of the paper is as fol- 
lows: Sect. 2 gives information on the selected sources and Sect. 3 
describes the data and analysis methods. Sect. 4 details the general 
spectral modelling procedure. The results are presented in Sect. 5. 
A general discussion is given in Sect. 6 while Sect. 7 offers the 
conclusions 



2 THE SOURCES 

2.1 The X-ray selected Radio Loud AGN sample from the 
Einstein Medium Sensitivity Survey 

Wolter & Celotti (2001) selected 39 Radio-Loud, broad line AGN 
from the X-ray EMSS sample (Gioia et al. 1990, Stocke et al. 
1991). Of these, 20 had a measured radio spectral index a,. < 0.7 
and their broad band properties were compared to those of classi- 
cal radio selected flat spectrum quasars. Their broad band indices 
(aro, a.ox) were found to be similar to those of FSRQ. This is at 
odds with the case of BL Lac objects, which occupy very different 
regions in the q^o, Q-ox plane depending on whether they derive 
from X-ray or radio selection. On the other hand, broad-band in- 
dices give only a global information, which may be ambiguous (see 
below). We therefore decided to exploit the unique capabilities of 
the Swift satellite to observe the ten X-ray brightest FSRQs of the 
EMSS sample in order to derive their X-ray spectrum together with 
simultaneous optical data allowing a reliable snapshot of the opti- 
cal to X-ray SED. Observations were performed within the filler 
program. The data gathered and analysis methods are described in 
Sect. 3. 



2.2 Controversial Sources 

In addition to the EMSS sources, we include in our analysis and dis- 
cussion four sources of controversial classifications put forward by 
other authors: two ultraluminous, high-redshift, broad line quasars 
recently discovered, SDSS J081009.94 + 384757.0 {z = 3.946, 
Giommi et al. 2007b) and MG3 J225155 + 2217(2 = 3.668, Bas- 
sani et al. 2007) have been claimed, on the basis of their broad-band 
spectral indices, to show a synchrotron peak in hard X-rays. If con- 
firmed, these claims would imply strong violations of the sequence 
scheme. 

Two other cases concern blazars of intermediate luminosi- 
ties, namely RX J1456.0 -|- 5048 {z = 0.478, Giommi 2008) and 
RGB J1629 + 401 (2 = 0.272, Padovani et al. 2002), which show 
synchrotron peaks in the X-ray band and, in addition, exhibit emis- 
sion lines, weak in the first source and pronounced, but narrow, in 
the second one. For these 4 objects, all the existing data obtained 
with Swift and, in the case of MG3 J225155 -|- 2217 also the data 
obtained with INTEGRAL were systematically reanalysed to pro- 
duce a data set as homogeneous and complete as possible. 



3 DATA ANALYSIS 

All the data presented here were analyzed in the same way to guar- 
antee a homogeneous treatment. The common procedures adopted 
for the analysis are described below. 
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3.1 Swift 

The data from all three instruments onboard Swift, namely BAT, 
XRT and UVOT (Gehrels et al. 2004), have been processed and 
analyzed with HEASof t v . 5.3.2 with the CALDB release of 
November 11, 2007. The observation log is reported in Tables [T] 

The X-ray Telescope XRT (0.2-10 keV, Burrows et al. 2005) 
data were analyzed using the xrtpipeline task, selecting single 
to quadruple pixel events (grades — 12) in the photon-counting 
mode. Individual spectra from multiple pointings, together with the 
proper response matrices, were integrated by using the addspec 
task of FTOOLS. The final count spectra were then rebinned in or- 
der to have at least 20 — 30 counts per energy bin, depending on 
the available statistics, and spectral fits were performed with simple 
power law models and galactic Nh- In only three cases additional 
absorption or a broken power law fit was required. The average 
spectral parameters are reported in Table |2] 

The sensitivity of the hard X-ray detector BAT (optimized for 
the 15 — 150 keV energy band, Barthelmy et al. 2005) is generally 
insufficient to detect extragalactic sources in short exposures, un- 
less for bright states. Therefore, in order to get meaningful signal- 
to-noise ratios we combined all the available pointings for a given 
source. All the shadowgrams from pointings of the same source 
were binned, cleaned from hot pixels and background (flat-field), 
and deconvolved. The intensity images were then integrated by us- 
ing the variance as weighting factor. No detection with sufficient 
signal-to-noise ratio has been found for any source with the avail- 
able exposures and upper limits at 3(j (corrected for systematic er- 
rors) in two energy bands (20 — 40 and 40 — 100 keV) have been 
calculated (see Table|3]l. 

Data from the optical/ultraviolet telescope UVOT (Roming et 
al. 2005) were analyzed using the uvotmaghist task with source 
regions of 5" for optical filters and 10" for the UV filters, while the 
background was extracted from a source-free annular region with 
inner radius equal to 6" (optical) or 12" (UV) and outer radius from 
30" to 60", depending on the presence of nearby contaminating 
sources. In the cases of MS 0402.0 - 3613, SDSS J081009.94 + 
384757.0, RX J1456.0 + 5048, and MG3 J225155 + 2217 the 
presence of a nearby source prevents the use of a background re- 
gion with annular shape; therefore, a circular source free region 
with 60" radius was selected. We added a 10% error in flux (corre- 
sponding to about 0.1 magnitudes) to take into account systematic 
effects. The summary of average magnitudes per filters is reported 
in TableH 



3.2 INTEGRAL 

The hard X-ray (20 - 100 keV) data on MG3 J225155 + 2217 
were obtained from the IBIS/ISGRI detector onboard INTEGRAL 
(Lebrun et al. 2003). They were analyzed with the Off-line 
Scientific Analysis (OSA) v. 7 . 0, whose algorithms 
for IBIS are described in Goldwurm et al. (2003), using the latest 
calibration files (v . 7.0.2). 

The source was observed serendipitously during the revolu- 
tions 316 (May 15 - 18, 2005) and 337 (July 17 - 20, 2005) 
for a total exposure of 388 ks. The results are reported in Table[3] 
The joint data from the INTEGRALIlSGm together with those of 
Swj/f/XRT, covering an energy range of 0.2 — 100 keV, can be still 
fit with a single power-law model with F = 1.41 ± 0.09 and ad- 
ditional absorption of TVg = {2.lt\\) x 10^^ cm"^, to mime the 
low-energy photon deficit. The fit gives a = 0.96 for 50 degrees 



of freedom. The results are consistent with the values for SwiftfKSI 
only, reported in Table[2] 



3.3 Two-point spectral indices 

Two point spectral indices aro, o-ox, ctrx between fixed 
radio/optical/X-ray frequencies have been computed for all 
sources, to facilitate comparison with other samples. 

The two-point optical/X-ray spectral index has been calcu- 
lated according to the formula (Ledden & O'Dell 1985): 



Ql2 = 



iog(g2/5i: 

log(!^2/;^l) 



(1) 



where 5*1 and S2 can be the radio (Vr ~ 5 GHz), optical (at 
Uo = 8.57 X lO^** Hz corresponding to 3501 A, U filter) and X-ray 
(at = 2.42 X 10^^ Hz corresponding to 1 keV) flux densities, 
respectively. The flux densities have been K-corrected by multi- 
plying for the factor (1 + z)"~^, where a is the spectral index of 
5*1, oc v'". An average spectral index a has been used for radio 
and optical frequencies, with values of ar — 0.5 and Oo = 1-38 
(Plan & Treves 1993). The and X-ray flux density have been 
measured from the spectral fitting of Swift/XKF data reported in 
Tabled 

The radio flux densities at 5 GHz have been extracted from 
neeQ or from radio catalogs available through HEASARCB When 
the 5 GHz flux density was not available (a very few cases), we have 
extrapolated the requested value from the 1.4 GHz flux density. 

The SwiftfUYOT flux densities have been calculated from the 
observed magnitudes, dereddened with the Av from NED and ex- 
trapolated to other frequencies by means of the extinction law by 
Cardelli et al. (1989), and converted by using the zeropoints given 
in the HEASof t package. 

The derived values are reported in Table |4] 



4 SED MODELS 

The observed SEDs of all sources were modelled following the 
theoretical scheme developed by Ghisellini, Celotti & Costamante 
(2002) and also used by Celotti and Ghisellini (2008, CG08), where 
a more detailed description can be found. 

Briefly, the observed radiation is postulated to originate in a 
single dissipation zone of the jet, described as a cylinder of cross 
sectional radius R and thickness (as seen in the comoving frame) 
AR' = R. We assume a conical jet with aperture angle ^pj — 0.1, 
the dissipation is assumed to always occur at i?diss ~ 10 R. 

The relativistic particles are assumed to be injected through- 
out the emitting volume for a finite time tj^j = AR' /c with total 
injected power in relativistic particles L(„j in the comoving frame 
of the emission region. The observed SED is obtained from the par- 
ticle distribution resulting from the cumulative effects of the injec- 
tion and radiative energy loss processes at the end of the injection, 
att — t'l^j, when the emitted luminosity is maximised. 

The finite injection duration, together with the magnetic and 
radiation energy densities seen by the particles, determine the en- 
ergy 7c above which the injected particle spectrum is affected by 
the radiative cooling processes. 

The particle injection function is assumed to extend from 7 = 



^ http://nedwww.ipac.caltech.edu/ 
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Figure 1. Spectral Energy Distributions (SED) of the 10 EMSS blazars in our sample, shown in order of increasing redshift. The models include the contribution 
of a putative accretion disk (dashed lines) described by a blackbody component peaking at ~ 10^^ Hz plus a power law component in X-rays produced by a 
hot corona. The non thermal jet emission (dotted lines) is due to synchrotron radiation (selfabsorbed at mm wavelengths) and IC scattering of photons in the 
broad hne region (see text). The solid lines represent the total emission. Light grey symbols (green in the electronic version) refer to archival data (taken from 
NED), darker symbols (red in the electronic version) represent the data analyzed in this paper The gray line indicates the expected flux sensitivity (Sir, 1 year 
of exposure) of the FGRT. 
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1 to 7max, with a broken power-law shape with slopes oc 7^ and 
oc 7"" below and above 7inj. 

If 7c < 7iiij (fast cooling regime) the resulting particle spec- 
trum iV(7) is given by 

A'"(7) oc 7~^'''''''"^; cx 7"^^; oc 7""*"; 

7 > 7inj; 7c<7<7inj; 7 < 7c (2) 

In the opposite case (slow cooling regime), radiative losses affect 
the spectral shape only at higher energies 7 > 7c > 7inj! yielding: 

Nipi) oc 7"^""''"'"'; oc 7~*; oc 7~'' 

7 > 7c; 7mj<7<7c; 7<7mj- (3) 

The relativistic electrons radiate via the synchrotron and the Inverse 
Compton mechanisms. The latter includes the SSC and the exter- 
nal Compton (EC) processes. The external radiation for the EC pro- 
cess is assumed to derive from the accretion disk emission scattered 
and reprocessed in the surrounding broad line region (Sikora et al. 
1994). 

With the distributions given above, the synchrotron and IC 
peaks in the SED are due to electrons with Lorentz factor, 7pcak, 
determined by the injected distribution and the radiative losses: in 
the fast cooling regime 7peak = 7inj as long as 7inj > 1. When 
7inj = 1, 7pcak = 7c (provided that s < 3). In the slow cooling 
regime, instead, 7peak = 7c (if s < 3). 

In most of the sources considered here the data indicate the 
presence of a prominent optical/UV excess that one can identify 
with the disk emission and then can be used to estimate Ldisk. The 
disk emission is assumed to be a black body, peaking at frequency 
I'ext ~ 10^^ Hz. Besides this black body UV component, we also 
consider the emission of the X-ray corona. This X-ray component 
is assumed to have ax = 1, to start a factor 30 below the UV bump, 
and to have an exponential cutoff at 150 keV. When no opticalAJV 
excess is apparent we still include a possible accretion disk (that 
can be considered as an upper limit). The external radiation field is 
derived assuming for the Broad Line Region a luminosity I/blr = 
O.ll/disk and a size -Rblr derived from the relation 

Rblr ~ 10"idLl45 cm, (4) 

where I/disk ~ 10*^idisk,45 erg s~^, in approximate agreement 
with Kaspi et al. (2007). Since [/blr = I/BLR/(47ri?|Lp^c), the 
radiation energy density of the BLR photons, as measured by an 
observer inside the BLR at rest with respect to the black hole, is 
constant, and equal to (7blr ~ 0.026 erg cm~"^. In the comoving 
frame, however, the radiation energy density will be enhanced by a 
factor. When -Rblr is smaller then the dissipation zone -Rdiss, 
we completely neglect any external radiation. 

The essential parameters of the model and their role can sum- 
marized as follows: 

• R, iinji s and 7inj determine the injected relativistic particle 
distribution; 

• the magnetic field B and the radiation energy density [/rad 
determine 7c, the energy of the cooling break, and thus the parti- 
cle spectrum after cooling. Umd includes the contribution of both 
the synchrotron photons and the external photons. In the comoving 
frame U^^^ ~ T^Ucxt, therefore also the bulk Lorentz factor F is 
involved. 

• the viewing angle 9, together with F, determines the Doppler 
(beaming) factor 5, and then the boosting of the entire SED; 



• the ratio F^f/cxt/B^ determines the EC to synchrotron lumi- 
nosity ratio. 

The SED peak frequencies (i^pcak.s, i^pcakjc) depend mainly 
on the shape of the particle spectrum. If 7c is high (low cooling) 
7pcak > 10* and !^pcak,s is above the optical band and the SSC 
dominates the high energy emission. If 7c is low (strong cooling) 
7pcak < 10-100 and then !^pcak,s is very much below the optical 
band, often below the self-absorption frequency, and the EC mech- 
anism produces the high energy component. 

Finally, the power carried by the jet in different forms (ra- 
diation, magnetic field, relativistic electrons, cold protons) can be 
computed from the model parameters as: 

Li ~ nR^V'^ficUi (5) 

where [/rad = f^r'ad ' Ub = B^/8tT, Uc = Ue < J > rTleC^, 

Up = npTJipC^ 

We refer to CG08 for a recent estimate and discussion of the 
jet powers for a large number of blazars. 



5 RESULTS 

5.1 The EMSS sample 

In all cases the X-ray spectra are well fit by simple power laws with 
photon indices ranging from 1.8 to 1.3 and the aox index, which in- 
dicates the relative intensity of the X-ray to optical emission, varies 
from 0.9 to 1.6. None of the objects is detected in the hard X-rays 
with BAT. 

In radio-quiet AGN, where the optical as well as X-ray emis- 
sion derive from an accretion disk, this ratio is on average aox — 
1.6 (see, e.g., Vignali et al. 2003). Thus, the lower values of aox 
for the EMSS sample, i.e. the high ratios of X-ray to optical fluxes 
(see Table |4} point to a significant contribution from the jet in the 
X-ray band. This conclusion is reinforced by the derived spectral 
indices in the X-ray band (ax — 0.3 — 0.7) which are unusu- 
ally hard for normal, even radio loud but not flat spectrum, quasars 
(ax — 0.8—1, see Grandi et al. 2006, Guainazzi et al. 2006). There- 
fore, the observed X-ray emission can be attributed to the combined 
contributions from the accretion disk corona plus a harder X-ray 
component deriving from Inverse Compton scattering of external 
photons (EC) off relativistic electrons in the jet. We recall that all 
the objects in the sample have Broad Line Regions providing seeds 
for the EC process. 

The optical and X-ray fluxes derived from Swift for the 10 X- 
ray selected radio-loud AGN are shown in Fig. [T] (in order of red- 
shift), together with historical radio and optical data, and theoretical 
models for the full SEDs. The high optical/UV fluxes and hard X- 
ray spectra indicate that the UV to X-ray SED have a concave shape 
as commonly observed for FSRQs. Therefore, in all cases the syn- 
chrotron peak frequency should fall below the UV range (red SED). 
On the other hand, the lack of data in the IR-submm region does not 
allow to set stronger observational constraints on its location. 

It is noteworthy that the two objects with the highest aox, 
MS 2141.2 + 1730 and MS 0521.7 + 7918, have relatively low 
radio luminosity (~ 10*^ erg s~^). This suggests that the jets in 
these sources could be either intrinsically weak or seen at inter- 
mediate angles, therefore not strongly beamed (see also Landt et 
al. 2008). Indeed, the unified model predicts that blazars at inter- 
mediate angle should exist (in large numbers). Sources with sublu- 
minous jets were not present in the samples used to construct the 
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Figure 2. Spectral Energy Distributions of MS 1340.7 + 2859, together 
with two possible fits. The one with the high energy emission above the 
GLAST detection curve (blue line), is the same as in Fig.[T] The parameters 
for the two fits, listed in Tab.|5] differ only slightly, but the predicted flux in 
the GLAST band is quite different. 



sequence, likely because the high threshold in radio flux selected 
only the most beamed sources. 

The X-ray selected AGN considered here, whose radio emis- 
sion was in some sense measured "a posteriori", have radio fluxes 
on average around 200 mjy and radio luminosities in the range 
one order of magnitude lower than the 
range covered by the FSRQs in the sequence. The optical fluxes, 
indicative of the accretion disk emission, correspond instead to lu- 
minosities of lO*^^*'^ erg s~^, well in to the typical quasar range. 
An even lower radio power interval was explored by Caccianiga 
and Marcha (2004), though without spectral information. 

MS 0828.7 -I- 6601 is extreme in the opposite way, in that the 
optical luminosity is lower than the X-ray one and the X-ray spec- 
trum is extremely hard (a^ = 0.3). The Galactic absorption col- 
irnin is A^H = 4.31 x 10^" cm~^ and the X~ray spectrum does not 
show any hint of additional absorption, suggesting that the source 
is intrinsically faint at optical wavelengths. 

The SED models shown in Fig. 1 (continuous/blue lines) were 
computed as described in Section 3. As noted above the main obser- 
vational constraints are that the synchrotron peak should fall below 
the UV range and that the jet contribution in X-rays should have 
a hard slope. Observational data below the optical and above the 
X-ray range are lacking so that the models cannot be more strongly 
constrained. According to our SED models the peak of the high en- 
ergy emission, always dominated by the EC process, falls around 
1-10 MeV, unusually low compared to known 7~ray blazars. 7- 
ray observations would be extremely important to confirm these 
models: the sensitivity curve of the Fermi GRT,for a 1 yr exposure, 
shown for comparison on each SED, indicates that several of these 
sources may be detected in the near future. 

The model parameters adopted are reported in Table |5] 

Since the lack of high energy (hard X-ray and 7-ray) data 
leaves some freedom in the choice of the fitting parameters we 
present also examples of alternative models for the SEDs. 
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Figures. SED of MS 1050.9-f 5413 andMS 1623.4-f 2712 together with 
two possible fits. For both sources the dashed line is the same fit as in Fig. 
[T] The solid line, instead, corresponds to assuming that the emission region 
is beyond the radius of the broad line region. Consequently, the EC process 
is completely neglected. This corresponds to less radiative cooling and a 
larger 7peak> and in turn this corresponds to a synchrotron peak at much 
larger frequencies. In this case, the high energy flux is completely due to 
the SSC process. 



In Fig.[2]two models for the SED of MS 1340.7 -f 2859, both 
fitting the optical and X-ray data equally well, are compared. The 
first (same as in Fig.[Tll predicts a 7-ray flux well above the Fermi- 
GRT sensitivity; the second predicts a lower 7-ray flux, just below 
the FermiGRT sensitivity limit. The parameters for the two models 
are quite similar (see Tab.[5]l. 

In Fig.[3]altemative SED models (sohd fines) for MS 1050.9+ 
5413 and MS 1623.4-f 2712 are shown. Both sources are character- 
ized by a relatively small accretion disk luminosity, hence a small 
^?BLR, a "normal" X-ray to optical luminosity ratio and a "nor- 
mal" X-ray slope (as opposed to e.g. MS 0828.7-1-6601, which has 
a relatively stronger, and harder. X-ray component). We hypothe- 
size here that the dissipation region of the jet is beyond the BLR 
(i.e. i?diss > i?BLR), contrary to what assumed for the models 
shown in Fig.[T] Then, in the absence of external photons, the cool- 
ing threshold 7c is high and the synchrotron emission peaks close 
to (but not above) the optical band, while the whole high energy 
emission from X-rays to 7-rays is produced by the SSC (instead of 
EC) process. Note that the injection spectrum is also different from 
the previous models, in that 7inj changes from 1 to 400 and 300 
for the two sources respectively (see Tab.[5]l. As a consequence the 
jet powers computed for the latter models are significantly reduced 
(Tab.llll. 

These examples illustrate the ambiguities in our models when 
Ux and ttox are "normal", and the accretion disk luminosity is rela- 
tively modest, suggesting a small broad line region. It is clear from 
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Fig.[3]that additional data, especially between the sub-mm and the 
optical region of the SED are needed to discriminate between these 
alternatives. In the following we discuss the results assuming in all 
the cases the parameters obtained with the standard external Comp- 
ton model. 

5.2 Controversial blazars 

5.2.1 Two high red-shift "red" blazars 

Two high redshift powerful blazars recently discovered 
(SDSS J081009.94 + 384757.0, Giommi et al. 2007b; 
MG3 J225155 + 2217, Bassani et al. 2007, Falco et al. 1998), 
both showing conspicuous broad lines, have been claimed to 
possibly exhibit a synchrotron peak in the X-ray band at variance 
with the sequence scheme. Our reanalysis of all the existing 
Swift and INTEGRAL data is summarized in the Tables and the 
results are shown in Fig. |4] For SDSS J081009.94 + 384757.0 
the Swift observations yield only upper limits in the optical UV 
range. Therefore the optical fluxes measured by the SDSS, though 
not simultaneous to the X-ray data, are also shown in Fig. |4] 
In both cases the hard X-ray spectra suggest a concavity of the 
optical to X-ray SEDs, pointing to an inverse-Compton origin 
of the X-ray/hard X-ray emission. Though both objects have 
aRx — 0.6, less than the conventional threshold value of 0.78 
taken to distinguish HBL from LBL, the "red" nature of their SEDs 
seems clear even with these limited data. 

The case of MG3 J225155 + 2217 discovered with INTE- 
GRAL is extraordinary: the X-ray to hard X-ray component domi- 
nates the optical emission by more than one order of magnitude. 
Attributing the latter component to synchrotron emission would 
require implausibly high values of both magnetic field and par- 
ticle energies. The SED of this source is very similar to that of 
Swift J0746.3 + 2548 (2 = 2.979) discovered by the BAT instru- 
ment on board Swift (Sambruna et al. 2006a). 

Clearly, a hard X-ray selection favors the discovery of objects 
with an extremely dominant hard X-ray component. It is neverthe- 
less interesting that blazars with such extreme dominance of the 
high energy emission component actually exist. 

Our model SEDs are shown in (Fig. |4ll For 
SDSS J081009.94 + 384757.0, although the wavelength 
coverage is poor, the hard X-ray spectrum measured by Swift 
strongly disfavors both a synchrotron or SSC interpretation of the 
X-ray emission. We therefore propose a model with synchrotron 
peak at low energies and a strong EC component. On the other 
hand the peak energy of the EC component is uncertain as in the 
case discussed above for MS 1340.7. 

For MG3 J225155 + 2217 the INTEGRAL data define well 
the lower energy side of the EC component: this branch of the SED 
reveals the low energy end of the electron distribution. The SED 
model is sufficiently determined by the data and the peak energy of 
the EC component is unlikely to be higher. 

5.2.2 Two luminous "blue" blazars 

The SEDs for RX J1456.0 + 5048 and RGB J1629 + 401 including 
all the data analysed here, are shown in Fig. [5] Here the data con- 
firm the "blue" nature of the SEDs as discussed already in Padovani 
et al. (2002) and Giommi (2008) respectively. In both cases the syn- 
chrotron peak falls in the UV to X-ray range as typical for "High 
energy peaked" BL Lacs (HBL). 

Indeed RX J1456.0 + 5048 can be classified as a BL Lac, 
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Figure 4. Spectral Energy Distributions (SED) of the 2 high-redshift 
blazars: MG3 J225155 + 2217 and SDSS J081009.94 + 384757.0. Same 
symbols as in Fig.[T] 

since its optical spectrum from the SDSS (see Fig.|6) shows a blue 
continuum and very weak emission lines, with equivalent width 
< 5 We note also that the radio source previously associated 
with this object is probably a misidentification: a bright, compact 
flat-spectrum radio source in the field, about 44 arcsec to the east 
of the SDSS BL Lac object, dominates the radio flux in this region. 
In the best image from the FIRST, at 1.4 GHz, the field source is 
168 miy, whereas the BL Lac object is only 9.5 mjy (T. Cheung, 
private communication). The latter value for the radio flux is shown 
in Fig. [5] RX J1456.0 + 5048 is highly luminous, but still compa- 
rable to the well known HBL PKS 2155 - 304 in the bright state 
recently observed (Foschini et al. 2007, 2008). 

The SED model for this object shown in Fig. |5] formally as- 
sumes i?diss > J?BLR, however there is no indication that a broad 
line region exists at all in this object. Also the possible contribution 
of an accretion disk is only an upper limit. The cooling threshold 7c 
is very high, as well as the break energy of the injection spectrum 
7inj which coincides with 

7pcak • 

The case of RGB J1629 + 401 is intriguing. Its optical spec- 
trum, also from the SDSS, is shown in Fig. |7] The emission lines 
are pronounced, but "narrow" {FW HM < 1500 km/s, Komossa 
et al. 2006), typical of Narrow-line Seyfert 1 galaxies. Estimates of 
the mass of the central BH are in the range 2 x lO'^ Mq (Komossa 
et al. 2006). This is unusual for a radio loud source and in particu- 
lar for blazars. However mass estimates for NLSyls may have to be 
revised (Marconi et al. 2008, Decarli et al. 2008). In fact the whole 
case of NLSyl and particularly that of radio loud NLSyl is the sub- 

According to the line measurements available in the SDSS server at 
|http://cas.sdss.org/astrodr6/en | 
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Figure 5. Spectral Energy Distributions (SED) of the 2 "blue" blazars with 
emission lines, RX J1456.0+5048 and RGB J1629+401. The plotted disk 
and corona spectra should be considered as an upper limit to the thermal 
emission in these blazars. Since the dissipation region is beyond the location 
of the BLR, their photons so not contribute to the inverse Compton spectra. 
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Figure 6. Optical spectrum of RX J1456.0 + 5048 from the Sloan Digital 
Sky Survey (SDSS). 

ject of active debate (e.g. Malizia et al. 2008, W. Yuan et al. 2008). 
Specifically RGB J1629 + 401 is only moderately radio loud and 
its X-ray emission spectrum showing a broken power law spectrum 
with spectral indices Ti = 3.4^°;^ and Ta = 2.21 ± 0.06 below 
and above 0.5 keV does not point unambiguously to a jet contri- 
bution. The X-ray emission could be driven by its NLSyl nature 
(Komossa et al. 2006) 

Nevertheless we modelled the SED of this source assuming 
a jet contribution and negligible radiative losses in the dissipation 



Figure 7. Optical spectrum of RGB J1629 + 401 from the Sloan Digital 
Sky Survey (SDSS). 

region (i?diss > i?BLR )- The radiation energy density is thus due 
to magnetic field and synchrotron photons only. The resulting SSC 
component has low luminosity. 

The parameters derived modeling the observed SEDs for all 
the objects are given in Table [S] The powers carried by the jets in 
various forms, derived from the models, are reported in Table[6] 



6 DISCUSSION 

6.1 Model parameters and correlations 

It is interesting to compare the (model independent) two point spec- 
tral index aox with quantities derived from the modelling for the 
EMSS sample, aox describes the ratio of the optical to X-ray emis- 
sion. For the "red" SEDs of this sample the optical is often domi- 
nated by the accretion disk emission while the X-ray contribution 
in excess of the accretion disk corona can derive from the jet, there- 
fore aox is a measure of their relative strength. 

We find from the modelling that aox correlates with the view- 
ing angle 6 (Fig. [8} as well as (inversely) with the Doppler factor 
5 (not shown). We recall that V and 5 enter the models in different 
ways so that both are determined in the fits: 6 is then derived from 
the first two. The two correlations are not independent as the value 
of r is almost constant ranging between 10 and 15. 

The aox -0 correlation agrees with the initial expectation that 
relatively weak jet contributions could derive from intermediate 
viewing angles. However aox also correlates inversely with the es- 
timated jet intrinsic power, so that the two effects, a lower degree of 
beaming or a lower intrinsic jet power, cannot be disentangled and 
probably coexist. On the other hand, aox does not correlate with 
Ldisk, indicating that it is the relative jet "strength", either due to 
the viewing angle or to the intrinsic power, that may drive these 
correlations. 

The most important parameter determining the shape of the 
SED is 'Ypcak, the energy of electrons radiating the peak syn- 
chrotron luminosity. The relation between 7peak and the radiation 
energy density as seen in the emission region frame, u^ad, is shown 
in Fig. [21 for the objects discussed here. Values derived for the 
blazars modelled in CG08 are also shown for comparison. The new 
objects fall well within the parameter sequence defined by the pre- 
vious sample. Therefore, despite the different selection criteria, X- 
ray selection as opposed to radio selection for FSRQ and selections 
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aimed at finding sources breaking the sequence trends, the spectral 
sequence holds in terms of physical parameters. 

The "red" objects tend to cluster at one end, with rather low 
7pcak- The clustering can result from the assumption (adapted from 
Kaspi et al. 2007) of an almost constant intrinsic radiation energy 
density in the BLR: since the magnetic field value is limited by 
the high Compton to synchrotron ratio, (/rad — UblbX^ (which 
varies by a factor 2 for V between 10 and 15) cannot be exceeded 
substantially. 

Blue SEDs are obtained for small (/rad, such that 7pcak = 7c 
is high, while the opposite is true for red SEDs. This is a result 
of the feed-back introduced by including radiative cooling in mod- 
elling the energy distribution of the relativistic electrons. However 
this is not the only important difference between the electron dis- 
tributions producing blue and red SEDs. Another major difference 
is at the lower electron energies: in fact, below 7pcak, the assumed 
electron distribution is extremely hard (oc 7^^) so that the total 
number of electrons is small and their average energy much higher 
than in the case of red SEDs (see Table 6). 



6.2 Jet powers 

The total jet power, Pjot , computed from the models (see Sect|4j 
is shown vs. the luminosity of the accretion disk (estimated from 
the data) in Fig. [TO] together with values previously derived for 
other blazars (Maraschi & Tavecchio 2003, Sambruna et al. 2006b). 
The two quantities tend to correlate for the red SEDs. Moreover 
the value of Pjet is of the order of 10 times the disk luminos- 
ity, corresponding to an approximate equality of the jet power and 
the accretion power (for an assumed radiative efficiency of the ac- 
cretion process of 10% appropriate for an optically thick flow). 
This represents on one hand an interesting consistency check on 
the models. In fact the bulk Lorentz factor, size, magnetic field, 
particle densities and spectra in the dissipation region are cho- 
sen to satisfy the SED constraints, including clearly the observed 
fluxes/luminosities, but the jet power is a global quantity that is 



Figure 9. The energy of electrons emitting at the peaks of the SED, 
7pcak a function of the energy density (radiative plus magnetic) as 
seen in the comoving frame. Oue EMSS sources are shown as empty cir- 
cles. The two squares correspond to the two "blue" FSRQs in our sample, 
while the two tilled triangles correspond to the two high redshift blazars 
SDSS J081009.94 + 384757.0 and MG3 J225155 + 2217. For compai'- 
ison, we show also all the blazars studied in Celotti & Ghisellini (2008) 
(grey asterisks). 



computed a posteriori. Moreover the new data and models increase 
significantly the number of sources in which a comparison of the 
jet power with the accretion power is possible and confirms the 
substantial balance between the two. 

The value of the power for jets with red SEDs depends quite 
strongly on the minimum energy of the injected electrons, 7min, 
which is often poorly constrained. For red SEDs the X-ray emis- 
sion derives from the EC process of relatively low energy elec- 
trons thus a lower limit can be obtained. Moreover in several cases, 
7min could be inferred from detailed fits of X-ray spectra showing 
a break in the X-ray range. Such breaks have sometimes been in- 
terpreted as absorption by cold gas (e.g. Cappi et al. 1997, Fiore 
et al. 1998, Fabian et al. 2001a,b, Bassett et al. 2004) but could 
represent instead the low energy "end" of the electron distribution 
(e.g. Tavecchio et al. 2007). With the latter interpretation 7miii is 
found to be in the range 1-10, thus supporting the high power val- 
ues (Tavecchio et al. 2000, Maraschi & Tavecchio 2003, CG08). 
These powers rest on the assumption of 1 cold proton per electron, 
but are independently required by the fact that the power directly 
radiated by the jet is larger than the power in relativistic leptons and 
magnetic field (CG08). 

For jets with blue SEDs the jet power derived from the present 
models is substantially lower, due to the smaller number and higher 
mean energy of particles in the assumed relativistic electron pop- 
ulation. The uncertainty on the amount of low energy electrons in 
the jet is however quite large because radiation from this branch of 
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Figure 10. The total jet power as a function of the disk luminosity, for 
the blazars studied in this paper. Filled red triangles are the high-2 blazars 
and filled blue circles are the "blue" blazars. For comparison, black open 
symbols report the data obtained for other blazars (Maraschi & Tavecchio 
2003, Sambruna et al 2006b). 



the electron distribution is not observable: its synchrotron emission 
is covered by the larger fluxes due to the outer regions of the jet and 
its SSC emission falls in the hard X-ray - soft gamma-ray region of 
the spectrum which is difficult to observe. Note that from Table 6 
the total kinetic power from the "blue" jets is hardly sufficient to 
account for the emitted radiation, thus the assumed spectra lead to 
a lower limit in the estimate of the jet powers. 

As a result, in the Pjct - Ldiak plane, the "blue" objects mod- 
elled here fall well below the Pjct = 10 Ldisk line. In the same 
diagram the grey open circles with upper limits to the disk lumi- 
nosities (horizontal arrows) represent BL Lacs from Maraschi & 
Tavecchio (2003) for which an electron distribution with slope 2 
extending down to a Lorentz factor of 1 was assumed below the 
peak energy. 

Despite the large uncertainties about the powers of the "blue" 
objects, it is interesting that the "red" blazars discussed here which 
result from independent selection criteria in the X-ray/hard X- 
ray band confirm the results previously obtained for radioselected 
FSRQ. 

Very different results were obtained by Nieppola et al. (2008) 
for a large sample of Blazars using beaming corrections derived 
from brightness temperatures at high radio frequencies using vari- 
ability to estimate the sizes of the emission regions. This technique 
gives information on the Doppler factors for the radioemitting re- 
gions in AGN. They find an anticorrelation between Doppler factor 
and peak frequency of the synchrotron component such that the ob- 
jects with higher peak frequency are less beamed. As a result the 
"intrinsic" powers of objects with large peak frequency (i.e. "blue" 
objects) are higher than those of "red" objects. 



While models of TeV emitting blazars require Doppler fac- 
tors as large as 10-50 (e.g. Rrawczynski et al. 2002, Konopelko et 
al. 2003, Finke et al. 2008, Aharonian et al. 2008, Tagliaferri et 
al. 2008), VLBI observations show very small, subluminal motions 
(e.g. Finer, Pant & Edwards 2008). These two evidences can be rec- 
onciled assuming that the jet strongly decelerates from the blazar 
region to the VLBI scale, from which most of the radio emission 
originates (Georganopoulos et al. 2003, Ghisellini et al. 2005). This 
would also fit in with the scenario in which FR I jets decelerate on 
moderate scale lengths while FR II jets may remain relativistic up 
to very large scales (e.g. Tavecchio 2007). Therefore the Doppler 
factors obtained by Nieppola et al. (2008) for the jet radio emit- 
ting regions are probably much smaller than the Doppler factors in 
the optical-to-7-ray emitting region especially in the case of high 
peaked BL Lac objects. 



7 CONCLUSIONS 

We have analysed and modelled new data for the SEDs of blazars 
including the first X-ray selected sample of Radio-loud Quasars and 
a few more blazars claimed to possibly challenge the "blazar spec- 
tral sequence". The first, model independent, conclusion is that the 
conventional separation of different "types" of blazars according 
to their two point spectral indices (Qro, ctoxyOirx) can be deceiv- 
ing, due to the complexity and variety of the SEDs: for instance the 
traditional criterion arx > (<)0.78 to distinguish between blue 
and red SEDs fails in the case of SDSS J081009.94 + 384757.0 
and MG3 J225155 + 2217 which have a^i 0.6, yet this is not 
due to synchrotron emission in the X-ray band but to an extremely 
dominant Compton component. Results from the Fermi Gamma- 
ray Telescope will be critical in confirming our interpretation 

The "blazar spectral sequence" concept still holds in terms of 
a parameter sequence. In Fig. 1 1 (left panel) the parameter 7pcak 
computed from SED models is plotted vs. total jet power also com- 
puted from the models. A "sequence" (i.e. correlation) is still ap- 
parent between these two quantities. We stress that 7pcak is related 
to the spectral shape while Pj^t is a global quantity related to the 
emitted luminosity. 

However (Fig. 11, right panel) 7pcak does not show a clear cor- 
relation with the observed radio luminosity, which was the basic 
quantity used in building the spectral sequence. This is likely due to 
the fact that at lower flux thresholds intrinsically luminous but less 
beamed objects enter the samples causing a mix of lower luminos- 
ity due to lower intrinsic power or lower beaming. This ambiguity 
could probably be solved with high resolution radio observations. 

The spectral sequence requires that different "types" of 
blazars, from FSRQ to low-energy peaked BL Lac Objects (LBLs) 
to HBLs, have SED peaks at increasing frequencies. This is ac- 
counted for by increasing values of 7pcak, derived in our models 
with the assumption that radiative cooling affects the electron dis- 
tribution becoming less important as the radiation energy density 
around the dissipation region decreases. However the particle dis- 
tributions needed to account for red and blue SEDs respectively 
differ not only at the high energy end but also at the low energy 
end. Blue SEDs require an electron population with higher "av- 
erage" energy, perhaps pointing to different injection/acceleration 
mechanisms. 

The present results show that even X-ray selection does not 
lead to find "blue" quasars as was the case for X-ray selected BL 
Lac samples (see also Landt et al. 2008). This strengthens the con- 
clusion that when a bright accretion disk is present the SED is 
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Figure 11. The energy of electrons emitting at the peaks of the SED, 7pcak ^s a function of total jet power (left, for comparison in gray data from CG08) and 
the radio luminosity at 5 GHz (right). Symbols as in Fig. 9. A "sequence" (i.e. correlation) is still apparent between 7poak and ^jot. while 7pcak does not 
show a clear correlation with observed radio luminosity. 



always "red". However, if the initial samples suggested luminos- 
ity/power as a fundamental parameter in determining the SEDs' 
systematics, the present study suggests that the fundamental con- 
dition that separates "red" SEDs from "blue" ones is the intensity 
of the radiation field in the emission region, which is dominated by 
reprocessed photons from an accretion disk, if present. The rela- 
tion with luminosity/power could be indirect, due to the fact that 
a bright accretion disk forms only for high values of the accre- 
tion rate in Eddington units (e.g. m > 10~^mcrit), while the ab- 
sence of it (RIAF, Radiatively Inefficient Accretion Flows) occurs 
for ra < lO^^rhcrit, thus at lower power /or a fixed mass of the 
accreting Black Hole. Thus, the original picture of blazar SEDs, as 
a one parameter family governed by luminosity, should be revised 
including the mass of the accreting Black Hole as an additional 
important parameter. A new scheme along these lines has been re- 
cently suggested by Ghisellini & Tavecchio (2008). 
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Table 1. Swift observation log. 
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Table 2. SwiftlXSI measurements. iVn (from Kalberla et al. 2005) and ATg are in units of 10^" cm~^; exposure is in kiloseconds; Ti indicates the photon 
index in the case of a single power-law model or the soft photon index in the case of a broken power-law model; for the latter, the columns r2 and Sbreak 
indicates the hard photon index and the break energy [keV]; the flux is in units of 10"^^ erg cm"^ 



Source 




An 


Exposure 






r2 




y-Zdof 


l-''> KJ keV 


MS 0152.4 + 0424 


1.132 


4.00 


10.0 




1.5 ±0.2 






0.58/8 


9.64 


MS 0232.5 - 0414 


1.439 


2.26 


15.4 




1.72 ±0.08 






0.93/25 


16.3 


MS 0402.0 - 3613 


1.417 


0.603 


27.7 


20 ±10 


1.76 ±0.09 






1.20/29 


18.3 


MS 0521.7 + 7918 


0.503 


6.93 


21.5 




1.9 ±0.2 






0.53/6 


3.70 


MS 0828.7 -h 6601 


0.329 


4.31 


9.3 




1.3 ±0.2 






0.12/5 


9.20 


MS 1050.9 + 5418 


0.995 


0.891 


11.2 




1.8 ±0.4 






0.86/4 


3.56 


MS 1234.9 + 6651 


0.852 


1.75 


32.8 




1.8 ±0.1 






1.37/19 


5.49 


MS 1340.7 + 2859 


0.905 


1.24 


12.1 




1.6 ±0.1 






0.42/12 


11.8 


MS 1623.4 + 2712 


0.526 


3.25 


15.6 




1.8 ±0.2 






0.63/10 


7.60 


MS 2141.2 + 1730 


0.211 


7.35 


43.0 




1.69 ± 0.06 






1.67/35 


15.9 


SDSS J081009.94 + 384757.0 
MG3 J225155 + 2217 


3.946 
3.668 


4.88 
4.90 


10.6 
23.1 




1.4 ±0.3 

1 41+0.10 






0.28/3 
1.05/42 


5.39 
31.0 



RX J1456.0 + 5048 0.478 1.60 4.8 2.2 ±0.1 0.77/25 17.7 

RGBJ1629 + 401 0.272 0.977 40.7 3.4t° J 2.21 ± 0.06 0.47±l°;g5 1.15/78 9.36 



Table 3. Svvi/f/BAT and INTEGRAL/lSGRl measurements. Exposures are in kiloseconds; fluxes are in units of 10^^" erg cm^^ s^^; upper limits are at 
3(7 level, after having taken into account systematic errors. The Crab spectrum taken as reference for conversion in physical units can be expressed with a 
power-law model with T = 2.1 and normaUzation at 1 keV equal to 9.7 ph cm~^ s~^ keV~^ (Toor & Seward 1974). 



Source 




Expositfe 


^20- 40 keV 


^40-100 keV 


MS 0152.4 + 


0424 


10.0 


< 2.1 


< 3.8 


MS 0232.5 - 


0414 


15.9 


< 2.8 


< 3.5 


MS 0402.0 - 


3613 


28.5 


< 1.7 


< 2.2 


MS 0521.7 + 


7918 


22.0 


< 2.4 


< 3.1 


MS 0828.7 + 


6601 


9.6 


< 3.2 


< 4.0 


MS 1050.9 + 


5418 


11.6 


< 2.6 


< 3.5 


MS 1234.9 + 


6651 


32.7 


< 1.5 


< 2.1 


MS 1340.7 + 


2859 


12.4 


< 2.7 


< 3.9 


MS 1623.4 + 


2712 


16.0 


< 2.5 


< 3.6 


MS 2141.2 + 


1730 


44.8 


< 1.8 


< 2.2 


SDSS J081009.94 + 384757.0 


11.0 


< 2.4 


< 3.8 


MG3 J22515E 


1 + 2217 


26.8 


< 2.3 


< 2.9 


MG3 J22515E 


1 + 2217* 


388 


0.12 + 0.03 


0.32 ± 0.07 


RXJ1456.0-f 


- 5048 


5.0 


< 3.3 


< 4.9 


RGB J1629 + 


401 


42.8 


< 1.3 


< 2.0 



* CVTEGiML observation. 
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Table 4. Swi/f/UVOT magnitudes (lower limits are at 3ij level) and the corresponding broad band spectral indices Qro. ckox and arx- These have been 
calculated by using the K-corrected fluxes at (3501 A), 1 keV and 5 GHz. See Sect. 3.3 for more details. 



Source 


V 


B 


U 


UVWl UVM2 


UVW2 




Q!ox 




MS 0152.4 + 0424 


18.6 + 0.2 


19.2 + 0.2 


18.5 + 0.1 






0.61 


1.05 


0.75 


MS 0232.5 - 0414 


16.6 + 0.1 


16.7 + 0.1 


15.5 + 0.1 






0.43 


1.29 


0.70 


MS 0402.0 - 3613 


17.5 + 0.1 


17.7 + 0.1 


16.6 + 0.1 






0.60 


1.01 


0.74 


MS 0521.7 + 7918 


17.1 + 0.1 


17.4 + 0.1 


16.6 + 0.1 






0.38 


1.48 


0.73 


MS 0828.7 + 6601 


> 19.7 


> 20.5 


> 20.3 






< 0.77 


< 0.91 


0.81 


MS 1050.9 + 5418 


19.3 + 0.3 


19.8 ± 0.3 


18.7 + 0.1 






0.57 


1.01 


0.71 


MS 1234.9 + 6651 


17.7 + 0.1 


18.2 + 0.1 


17.2 + 0.1 






0.48 


1.22 


0.71 


MS 1340.7 + 2859 


16.7 + 0.1 


17.1 + 0.1 


16.0 + 0.1 






0.39 


1.34 


0.69 


MS 1623.4 + 2712 


18.8 + 0.2 


18.9 + 0.1 


18.2 + 0.2 






0.60 


1.05 


0.75 


MS 2141.2 + 1730 


16.2 + 0.1 


16.4 + 0.1 


15.2 + 0.1 


15.1 + 0.1 


15.0 + 0.1 


0.49 


1.56 


0.83 


SDSS J081009.94 + 384757.0 


> 18.5 


> 19.4 


> 19.2 


> 18.7 > 18.9 


> 19.5 


< 0.34 


< 1.03 


0.56 


MG3 J225155 + 2217 






19.3 + 0.1 


18.9 + 0.2 18.4 + 0.1 


18.5 + 0.1 


0.52 


0.74 


0.59 



RX J1456.0 + 5048 
RGB J1629 + 401 



18.6 + 0.3 19.2 + 0.2 18.4 + 0.2 17.9 + 0.2 17.9 + 0.2 17.9 + 0.1 0.42 0.69 0.55 
17.9 + 0.2 18.2 + 0.2 17.2 + 0.1 16.8 + 0.1 16.5 + 0.1 16.5 + 0.1 0.35 1.26 0.66 
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Table 5. The input parameters of the model for oiffblazars. (1) Source name; (2) redshift; (3) radius R of emitting region [10^^ cm]; (4) intrinsic injected power 
[10''^ erg/s]; (5) bulk Lorentz factor; (6) viewing angle; (7) magnetic field intensity [Gauss]; (8) minimum random Lorentz factor of the injected particles; 
(9) maximum random Lorentz factor of the injected particles; (10) Tpeak! (H) spectral slope of particles above the cooling break; (12) disk luminosity 
[Iff*^ erg/s]; (13) radius of the BLR [lO^'' cm]. (14) random Lorentz factor of the electrons cooling in AR'/c. 



Source 


2 


R 




r 





B 


Tinj 


Tmax 


Tpeak 


n 


id 


^BLR 


7c 


Note 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


U) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


MS 0152.4+0424 


1.132 


10 


1.5e-2 


14 


3.3 


2.4 


1 


3e3 


18 


3.2 


5 


220 


11 




MS 0232.5-0414 


1.439 


18 


3.0e-2 


13 


4.0 


7 


1 


6e3 


50 


3.3 


100 


l.e3 


6 




U4U2.0-3D13 


1.41 / 




O la o 

1.5Q-1 


1 1 


4.U 


o 




A 

4 


ie4 


1 1 


3.4 


Ar\ 
4U 


/OO 


1 1 




MS 0521.7+7918 


0.503 


10 


l.Oe-2 


13 


6.0 


2 


150 


8e3 


150 


3.5 


7 


300 


18 




MS 0828.7+6601 


0.329 


2 


6.0e-5 


14 


3.2 


3 


1 


le4 


37 


3.45 


0.1 


25 


37 




MS 1050.9+5418 


0.995 


10 


1.5e-3 


12 


3.5 


5 


1 


le4 


21 


3.2 


2 


180 


21 


Rdiss < RbLR 






18 


2.5e-3 


12 


3.5 


0.12 


400 


2e5 


7e3 


3.5 


2 


140 


7e3 


^diss > ^?BLR 


MS 1234.9+6651 


0.852 


9 


7.0e-3 


12 


3.5 


3 


15 


5e3 


25 


3.6 


6.5 


320 


25 




MS 1340.7+2859 


0.905 


14 


2.5e-2 


13 


4.0 


3 


70 


3e3 


70 


3.6 


29 


500 


8 


L~, below GLAST 






14 


2.0e-2 


14 


4.5 


5 


40 


3e3 


40 


3.7 


25 


550 


9 


above GLAST 


MS 1623.4+2712 


0.526 


10 


1.6e-3 


11 


5.0 


3 


1 


le4 


37 


3.3 


1 


150 


37 


^diss < ^BLR 






15 


3.6e-4 


11 


3.5 


0.1 


300 


le4 


2.1e4 


3.6 


1 


100 


2e4 


^diss > ^BLR 


MS 2141.2+1739 


0.211 


10 


9.0e-4 


12 


5.0 


6.5 


20 


1.5e4 


20 


2.4 


4.5 


200 


12 




SDSSJ081009.94+384757.0 


3.946 


16 


0.1 


14 


3.0 


5 


25 


1.5e3 


25 


3.8 


100 


1000 


6 




MG3 J225155+2217 


3.668 


20 


0.2 


15 


2.5 


0.7 


35 


le3 


35 


3.8 


21 


500 


6 




RX J1456.0+5048 


0.478 


15 


2.0e-4 


15 


3.0 


0.3 


3.5e4 


5e5 


3.5e4 


3.4 


0.2 


45 


le4 


^diss > SbLR 


RGB J1629+401 


0.272 


10 


l.le-4 


10 


4.0 


1.5 


7e3 


le5 


7e3 


3.5 


0.25 


50 


930 


^diss > ^BLR 



Table 6. Kinetic powers and Poynting fluxes (all in units of 10 erg/s) — (1) Source name; (2) Total (synchrotron + IC) radiative power Lr; (3) Synchrotron 
radiative power Ls\ (4) Poynting flux Lb ; (5) Kinetic power in emitting electrons Z/e; (6) Kinetic power in protons Lp, assuming one proton per electron; (7) 
Average random electron Lorentz factor (7). 



Source 




U 


Lb 






(7> 


Note 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


MS 0152.4 + 0424 


3.15 


0.04 


0.42 


1.11 


253 


8.1 




MS 0232.5 - 0414 


4.92 


0.73 


10.03 


0.68 


152 


8.2 




MS 0402.0 - 3613 


2.29 


0.33 


4.16 


2.90 


1301 


4.1 




MS 0521.7 + 7918 


1.73 


0.04 


0.25 


0.29 


31 


17.6 




MS 0828.7 + 6601 


1.20e-2 


1.53e-4 


2.64e-2 


3.24e-2 


21 


2.9 




MS 1050.9 + 5418 


0.20 


2.67e-2 


1.25 


0.24 


129 


3.4 


-Rdiss < L^BLR 




0.39 


2.90e-2 


2.51e-3 


1.27 


14.4 


162 


^diss > ^BLR 


MS 1234.9 + 6651 


1.02 


2.71e-2 


0.39 


1.10 


236 


8.6 




MS 1340.7 + 2859 


4.32 


9.73e-2 


1.11 


0.78 


142 


10 


below GLAST 




3.80 


0.195 


3.59 


1.22 


254 


8.8 


L-,1 above GLAST 


MS 1623.4 + 2712 


0.19 


1.50e-2 


0.41 


3.66 


203 


3.3 


i?diss < RbLR 




5.60e-2 


6.0e-3 


l.Oe-3 


0.96 


15 


117 


^diss > ^BLR 


MS 2141.2 + 1739 


0.12 


1.38e-2 


2.27 


5.52e-2 


12.4 


8.2 




SDSS J081009.94 + 384757.0 


19.79 


0.44 


4.69 


7.10 


1996 


6.5 




MG3 J225155 + 2217 


50.05 


2.12e-2 


0.16 


12.77 


3320 


7 




RX J1456.0 + 5048 


3.26e-2 


2.13e-2 


1.70e-2 


1.17e-2 


7.1e-3 


3049 


^diss > ^BLR 


RGB J1629 + 401 


1.13e-2 


1.03e-2 


8.39e-2 


2.23e-3 


9.3e-3 


442 


^diss > Z^BLR 



